UCN upscattering rates in a molecular deuterium crystal 
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A calculation of ultra-cold neutron (UCN) upscattering rates in molecular deuterium solids has 
been carried out, taking into account intra-molecular excitations and phonons. The different molec- 
ular species ortho-Z)2 (with even rotational quantum number J) and para-D2 (with odd J) exhibit 
significantly different UCN-phonon annihilation cross-sections. Para- to ortho-D2 conversion, fur- 
thermore, couples UCN to an energy bath of excited rotational states without mediating phonons. 
This anomalous upscattering mechanism restricts the UCN lifetime to 4.6 msec in a normal-D2 solid 
with 33% para content. 
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The low density of ultra-cold neutrons (UCN) avail- 
able using conventional cold moderators at nuclear re- 
actors has long been the main constraint in the pursuit 
of high precision measurements of neutron /3-decay with 
UCN 1^. To pursue the possibility of utilizing this sim- 
ple hadronic system for tests of weak interaction theo- 
ries, Golub and Pendlebury ||^ proposed a way to in- 
crease UCN production through the exchange of energy 
between a cold neutron bath and the phonons in cer- 
tain cold moderators such as super-fluid "^He, with large 
neutron scattering cross-sections and small neutron ab- 
sorption cross-sections. Superthermal UCN sources ex- 
ploit the characteristics of low temperature substances, 
in which a large number of phonon modes are available 
for neutron downscattering, while the number of phonons 
present which can upscatter UCN is suppressed. Ideally, 
the density of UCN produced in such a source is limited 
either by nuclear absorption in the moderator as is the 
case for solid deuterium, or the neutron lifetime itself as 
is the case for superfluid "^Hc. 

A solid D2 UCN superthermal source Q is under devel- 
opment at the Los Alamos Neutron Science Center. Pre- 
liminary estimates jij promise gains in the available UCN 
density as high as two orders of magnitude over existing 
UCN facilities. To characterize its performance, calcula- 
tions of essential physical parameters, such as scattering 
cross-sections, UCN residence time, etc., with more de- 
tailed and accurate models than those presently available 
are in increasing demand. In this report, we demonstrate 
that in the presence of even very small concentrations of 
para-D2 (with total nuclear spin 1) can dominate the 
UCN upscattering rate, overwhelming the usual phonon 
annihilation mechanism. This results in greatly reduced 
UCN lifetimes in the solid and orders of magnitude re- 
ductions in the achievable UCN density. 

The deuterium molecule is a two-body system with 
the quantum properties of identical bosons. Its nuclear 
spin wave-function couples to molecular rotational states 
with same parity to preserve the symmetry of the wave- 



fmiction under permutation of identical particles. Koppel 
and Young calculated the neutron scattering cross- 
section of this molecular system, taking into consider- 
ation induced transitions between molecular rotational 
and vibrational states. In their formulation, an inco- 
herent approximation was used and translational coor- 
dinates were assumed to commute with inter-molecular 
degrees of freedom. The derived double differential cross- 
section in which the interference term is neglected has 
the form, for UCN with incident wavenumber k and final 
wavenumber fc', 
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where fin is the momentum transfer of the scattered neu- 
tron, Pjs is the population of the initial molecular state 
with a total nuclear spin S and rotational quantum num- 
ber J, Ej = 7meV x J(J -I- l)/2 is the rotational spec- 
trum, fiuj is the inter-molecular vibrational energy with n 
characterizing the number of vibrational energy quanta, 
and C(JJ'Z; 00) is a Clebsh-Gordon coefficient. Ani is 
defined as an integral over the orientation of a molecule, 
i.e.. 
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with a = O.74A the equilibrium separation distance of 
the D-D bond, /i the cosine of the inclination angle of the 
molecular axis from the z axis of a reference Euclidean 
coordinate system, and Pi the Legendre polynomial of 
order I. 
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TABLE I. Intrinsic scattering cross-sections Sjji associ- 
ated with different rotational transitions [6]. 



Sjj' 


Even-J(Ortho) 


Odd-J(Para) 


Even-J 
Odd-j' 


aloh + faLc = 6.687/47r 
faLc = 0.765/47r 


f aLc = 1.530/47r 
ciLh + 3a?nc = 6.f02/47r 



The input parameter Sjj' (in units of barns) in (||) for 
transitions between different rotational states is deduced 
and listed in table |. Note here that only the incoherent 
scattering length amc of a bound nuclide contributes to 
the ortho(even J)/para(odd J) conversion. 

In a D2 solid, the populations of the even-J(ortho 
state) and the odd- J(para state) are typically determined 
by the ortho/para population of the gas phase before the 
D2 is frozen into the solid. The self-conversion between 
these two species in the solid phase into a thermal Boltz- 
man distribution is extremely slow compared with the 
time scale of experiment. For example, a room temper- 
ature equilibrium D2 is a mixture of 67% ortho-D2 and 
33% para-D2. The conversion rate in the solid is mea- 
sured to be 0.06%/hr§, requiring about 7 months to re- 
duce the para content to 1.65% from 33%. On the other 
hand, relaxation to thermal distributions is rapid among 
ortho- and para- species themselves. Consequently, in 
low temperature circumstances relevant to super-thermal 
solid deuterium source(T<20K), only ground states( J — 
for ortho; J = 1 for para) are present. 

In the case of neutrons scattered off a low tempera- 
ture crystal with well-defined lattice structures, harmonic 
solid correlation functions should be applied to the trans- 
lational part of (Jl ) . Following the standard treatment of 
lattice dynamics we perform a phonon expansion 

< exp{-iK ■ ri}eXp{iK ■ r';/(t)} >Trans 

— exp{—2W{K)}exp{< k ■ uik ■ ui'{t) >}, 
= exp{~2W{Hi)}[l+ < K ■ uiK ■ ui'{t) > 

+ 0{< K ■ UlK ■ Ul,{t) >^)]. 

An overall Debye- Waller factor is extracted in front, and 
only the first two terms are left for discussion, yield- 
ing the zero and one phonon exchange processes, respec- 
tively. 

Unlike conventional applications of elastic solid cor- 
relation functions, the first term in the phonon expan- 
sion(zero phonon term) coupled to molecular internal 
energy states not only gives rise to UCN energy tran- 
sition(mainly upscattering), but overwhelms any phonon 
contributions when para-D2 is present. The conversion of 
para- into ortho- molecules provides direct energy trans- 
fer to UCN. This scattering cross-section without phonon 
couplings has the simple form: 
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where is a spherical Bessel function of order one 

The increase of the neutron momentum is definite, 

i.e., 



k' = y/2mnAEiQ/h. 
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A conversion energy Ai?io of 7 meV gives k' a value 
of 1.84x 10^°m~^. The momentum transfer k can be 
well-approximated by k' for UCN (fc' >> kucn—i-'27x 
10^m~^), and the above differential cross-section is 
isotropic, making the integration of (^) straightforward. 
A Debye- Waller factor originating from the uncertainty 
of positions of lattice sites, reduces the amplitude by a 
factor of 0.76. The temperature independently] total cross 
section of UCN upscattering ctio is calculated to be 31 
barns. This is at least of an order of magnitude larger 
than the phonon annihilation cross-section in a 4K solid. 
The rate of loss of UCN in the solid is 



Pucn = Wfi/V, 

— Pucn ucnf-' \ 
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Here p' is the density of para-D2, taken to be 3x10^^ 
cm^''. The corresponding upscattering time T„p of UCN 
in a pure para-D2 molecular solid is therefore 
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= 1.5 msec. 



For normal-D2 which retains the room temperature equi- 
librium ortho/para ratio, the upscattering time due to 
the spin relaxation of para species is 4.6 msec! 

To estimate the phonon upscattering rates, we approx- 
imate the SD2 hcp/fcc lattice as a cubic lattice to simplify 
the treatment of polarization anisotropics. The expres- 
sion of the incoherent double differential cross-section in- 
volving one phonon exchange is 
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where the energy of phonon iSp/i = e+AE{J' — > J)+nfiuj, 
complying with the law of conservation of energy. Pos- 
itive and negative values of Eph correspond to single 



*The only temperature dependence comes from a roughly 
5% decrease in the Debye- Waller factor between 4 and 18K. 
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FIG. 1. UCN upscattering cross-section vs. temperature 
of solid D2. The one-phonon annihilation cross-section in an 
ortho-D2 solid (solid curve) and in a para-D2 solid(dashed 
curve) are plotted. The dashed-dagger line is the temper- 
ature independent UCN upscattering cross-section involving 
J=l-^0 relaxation not coupled to phonons in a para-D2 solid. 



phonon creation and annihilation, respectively. Z(E) 
represents the normahzed phonon density of states, and 
n{E) the occupation number of phonons with energy E. 
With a simple Debye model, in which 
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and the Debye temperature Te(110K for D2) is the only 
parameter, a double integration of (^) with initial en- 
ergy of UCN (see Figure |l|) reproduces the upscattering 
cross-sections in ortho-deuterium calculated by Yu, Ma- 
lik and Golub in whose treatment rotational tran- 
sitions were not considered. In an ortho-D2 solid, the 
J = — > process dominates the upscattering. Even 
though the J = ^ 1 transition is energetically al- 
lowed through coupling to a phonon, the cross-section is 
kinematically suppressed by the smaller final phase space 
of upscattered neutrons. Para-D2 has a distinguishably 
larger one-phonon annihilation cross-section than the or- 
tho species. The origin of difference is again related to 
the J = 1 ^ relaxation channel. This provides UCN 
with additional energy to scatter into a larger volume of 
phase space, and secondly, it couples UCN to high en- 
ergy phonons with large density of states, and is thus 
less restricted by the availability of phonon modes than 
the J = 1 ^ 1 process. However, it is still suppressed by 
its small coupling to the phonon field, as opposed to the 
zero phonon term. 

In summary, para-deuterium has a spin relaxation 
channel in which its conversion energy of 7 meV can be 
released to UCN, resulting in a temperature-indepedent 
short UCN lifetime of 4.6 msec in a normal-D2 solid. 



Elimination of the para-D2 is necessary to achieve UCN 
lifetimes comparable to the nuclear absorption time in 
solid D2. 

The authors wish to acknowledge the generous support 
of the National Science Foundation through grant NSF- 
9807133 and the Department of Energy through LDRD- 
8K23-XAKJ funds. 



[1] 
[2] 

[3] 

[4] 

[5] 

[6] 

[7] 
[8] 
[9] 



P.R. Huffman et al, Nature, 403, 62(2000). 

R. Golub and J.M. Pendlebury, Phys. Lett., 53A, 

133(1975). 

R. Golub, D.J. Richardson, and S.K. Lamoreaux, Ultra- 
cold neutrons (Adam Hilger, Bistol, 1991), 66. 
Z-Ch. Yu, S.S. Malik and R.Golub, Z. Phys. B - Condensed 
Matter, 62, 137-142(1986) 

J. A. Young and J. Koppel, Phys. Rev., 135, 3A, 
A603(1964). 

These scattering lengths reproduce the calculation pre- 
sented by R.E. MacFarlane, Los Alamos National Labora- 
tory Report No LA-12639-Ms (ENDF 356), 1994. 
I.F. Silvera, Rev. Mod. Phys., 52, No. 2, Part I, 393- 
452(1980). 

S.W. Lovesey, Theory of neutron scattering from con- 
densed matter (Clarendon press, Oxford, 1984), Vol I, 98. 
N.M. Temme, Special functions - An introduction to the 
classical functions of mathematial physics (John Wiley & 
Sons, New York, 1996), 399. 



3 



